Static and dynamic shear viscosity of a single layer complex plasma 
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The static and dynamic (complex) shear viscosity of a single layer dusty plasma is measured by 
applying, respectively, a stationary and a periodically modulated shear stress, induced by the light 
pressure of manipulating laser beams. Under static conditions we observe a decrease of the viscosity 
with increasing shear rate, the so-called shear-thinning behavior. Under oscillating shear both 
the magnitude and the ratio of the dissipative and elastic contributions to the complex viscosity 
show strong frequency dependence, as the system changes from viscous to elastic in nature with 
increasing excitation frequency. Accompanying molecular dynamics simulations explain and support 
the experimental observations. 



I. INTRODUCTION 

Transport properties of strongly coupled complex (or 
dusty) plasmas, have attracted considerable attention 
during the past decade (see e.g. [l|). The interaction 
of the dust particles in these systems may be well ap- 
proximated by the Debye-Hiickel (or Yukawa) model po- 
tential $(r) = Q exp[—r / Xd] / Aireor, where Q is the 
dust particle charge, is the Debye screening length. 
Equilibrium Yukawa systems can be fully characterized 
by two dimensionless quantities: the Coulomb coupling 
parameter T = Q"^ / [AireQak-QT), where T is the temper- 
ature, and the screening parameter k = a/Xo- 

Viscosity, the measure of the plastic response of matter 
(primarily liquids and soft matter) to applied forces, is 
a central quantity in rheology. For a Newtonian fluid a 
constant shear viscosity r/ relates the shear stress a to the 
shear rate 7 = dvx/dy (velocity gradient) as cr = —777 — 
-rj{dvx/dy). 

Continuum hydrodynamics successfully uses the con- 
cept of viscosity, usually as input parameter, as an in- 
trinsic property of the material under investigation, in 
the Navier-Stokes equation. One has to note, however, 
that the Newtonian concept of viscosity is applicable only 
(i) at small shear rates, (ii) long length scales, and (iii) 
at low frequencies. In many physical systems, these con- 
ditions are clearly violated 0|- 

Regarding complex plasmas, the first experiment on a 
single layer (2D) dust cloud with control over the applied 
shear, via adjusting the power of a manipulating laser 
beam, was carried out by the group of Lin I in 2001 Q. 
Similarly, an experiment makiiig use of a single shear- 
ing laser beam was reported in [J] . In these experiments 
a sheared velocity profile was created around the beam. 
In the experiment reported in Q two displaced, parallel 
counter-propagating laser beams were used to realize a 
planar Couette configuration in a 2D dusty plasma layer. 
In another experiment the non-Newtonian behavior of a 
3D complex plasma in the liquid state was identified by 
Ivlev et al. [y]. Motivated by these pioneering studies. 



more recently, detailed dusty plasma experiments demon- 
strated the presence of viscoelastic response and re- 
vealed the wave-number dependence of the viscosity in 
2D 8]. Experiments in the crystalline phase have iden- 
tified the slipping of individual crystal lines to be the 
gimary mechanism for relaxing an applied shear stress 

Complementing the experimental efforts, the self- 
diffusion shear viscosity [T^-flsj and thermal 
conductivity coefficients [l^, [l^ have been derived in a 
number of simulation studies, both for three-dimensional 
(3D) and two-dimensional (2D) settings. In [g^], besides 
calculations of the "equilibrium" (7 — 0) static viscosity, 
predictions for the shear-thinning effect (typical for com- 
plex molecular liquids) were given at high shear rates. 
The frequency dependence of the complex shear viscos- 
ity, which combines the dissipative and the elastic com- 
ponents of the complex response of matter to oscillating 
shear stress as r]{ui) = ri'(uj) — ir]"{uj), was computed for 
3D Yukawa liquids in [21| . Fundamental questions about 
the existence or nonexistence of well defined transport 
coefficients in 2D were addressed in |22i] . 

Here we present laboratory dusty plasma experiments 
in which either a static or a periodic shear is applied on a 
single layer dusty plasma in the strongly coupled regime. 
During the evaluation of the experimental data we adopt 
an earlier experimental method @ as well as specific 
methods used so far only in molecular dynamics (MD) 
simulations. Our experimental results are supported by, 
and are combined with simulations and theoretical calcu- 
lations. To obtain the most complete information about 
the complex plasma single layer we perform three differ- 
ent experiments on the same dust cloud: (1) analysis of 
the thermally excited waves, without any applied shear, 
to determine the principal system parameters, (2) ap- 
plying a static shear (cjgh — 0) to investigate stationary 
flows at high shear rates, and (3) applying a periodic 
shear (wgh > 0) to obtain the frequency-dependent com- 
plex viscosity. 

The paper is structured as follows. Section II describes 
the details of the experimental apparatus. Section III 



discusses our experiments and the data evaluation, as 
well as the connection with simulation data. Section IV 
summarizes the work. 



II. EXPERIMENTAL APPARATUS 

Our dusty plasma experiments are carried out in a cus- 
tom designed vacuum chamber with an inner diameter of 
25 cm and a height of 18 cm. The lower, powered, 18 cm 
diameter, flat, horizontal, stainless steel electrode faces 
the upper, ring shaped, grounded aluminum electrode, 
which has an inner diameter of 15 cm and is positioned 
at a height of 13 cm. 

The experiments are performed in an argon gas dis- 
charge at a pressure p = 1.2 ± 0.05 Pa, and at a steady 
gas flow of ~ 0.01 seem. We apply 13.56 MHz radio 
frequency excitation of ~ 7 W power, to establish the 
discharge in which the melamine-formaldehyde micro- 
spheres with a diameter d = 4.38 ± 0.06 fim and a mass 
m = 6.64 X 10~^^ kg are levitated. For illumination of 
the particle layer we use a 200 mW, 532 nm laser, the 
light of which is expanded and enters the chamber from 
a side window. Our CCD camera has a resolution of 1.4 
Megapixels and captures snapshots of the particle en- 
semble through the top window of the chamber at 29.54 
frames per second acquisition rate. The camera is sensi- 
tive only at the illuminating laser wavelength, due to the 
application of an interference filter centered at 532 nm 
wavelength. 

The average dust particle number in the field of view 
is ^ 2500, while the total particle number in the dust 
cloud is about ~ 15000. During the evaluation of the 
raw images identification and position measurement of 
the particles is performed using the method described in 

M- 

The scheme of the optical setup used in the measure- 
ments to apply an external shear on the dust layer, is 
shown in Fig. [TJ The light of a 200 mW red diode laser 
(L) is split into two parallel beams of equal intensity and 
perpendicular linear polarization. One of these beams is 
directed towards another polarizer (P2), while the beam 
with the perpendicular polarization reaches P2 after be- 
ing reflected from a mirror. By rotating P2 by means 
of a motor the intensities of these two beams, passing 
through P2, can be modulated harmonically, with 180° 
phase shift relative to each other. These beams are fur- 
ther guided by additional mirrors, which are aligned to 
ensure that the two beams propagate inside the chamber 
in opposite directions, but on a common axis. This axis 
is adjusted to be horizontal and its vertical position is 
set in a way that the beams, which are focused to have a 
diameter of about 0.5 mm, lie within the dust layer. 

This setup makes it possible to exert a periodic force 
on the particles, which, in the center of the beam(s) has 
a form Fx{t) = Fq sin(wshi), with the direction x defined 
to coalesce with the beams, as shown in Fig. [T] The 
amplitude Fq can be tuned by the power of the laser 
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FIG. 1. (color online) Scheme of the optical setup used to 
generate the alternating shearing laser beams. Main parts: L 
- 200 mW 650 nm diode laser; PI - polarizing beamsplitter; 
M - mirror; P2 - linear polar filter mounted on a ball bearing 
and rotated by a motor. Beams with different arrow styles are 
modulated antiphase. In the reference beam line: R - 10 mW 
532 nm diode laser, M - mirror. 



(L), while the frequency Wsh = 2t:/T is defined by the 
half-period of rotation (T) of P2. As the recording of 
the images by the camera is not synchronized with the 
rotation of the polarizer, the angular position of the latter 
has to be measured. This is accomplished by a reference 
beam of a green diode laser (R). The polarized beam 
of this laser passes through P2 as well, and when P2 is 
rotated, the beam intensity is modulated. The beam is 
directed onto the surface of the lower electrode of the 
discharge, and, as the wavelength is the same as that of 
the laser illuminating the particle layer, the beam spot 
is picked up by the camera on the images of the particle 
layer. Measuring the intensity of this reference beam in 
the recorded images allows an accurate determination of 
the angular position of P2. 

In the experiments aimed at the study of the static 
viscosity P2 is not rotated, it is aligned to a position, 
which maximizes the intensity of one of the beams. 



III. EXPERIMENTS 

This section is divided into three parts. Section III. A 
explains the "calibration" of the system, which was aimed 
at the determination of the basic characteristics of the 
dust layer (particle charge, screening parameter of the 
potential, plasma frequency), needed in the subsequent 
data analysis. Section III.B discusses the measurements 
of the static viscosity covering a wide domain of applied 
shear rates. These measurements quantify shear thinning 
at elevated shear rates. A brief introduction of the simu- 
lation method is presented in this section. Section III.C 
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describes the experiments carried out with a harmonic 
perturbation, which reveal the viscoelastic behavior of 
the investigated system. 



A. Calibration 

The aim of this experiment is the determination of the 
basic parameters of our particle suspension, which will 
be crucial in the forthcoming experiments aimed at the 
measurement of viscous properties. Having calibrated 
the optical system the Wigner-Seitz radius is determined 
from the average particle separation. We find an areal 
density n — 3.21 mm~^ and the Wigner-Seitz radius to 
be a = \/ ^/wri = 0.315 mm. Fitting the particle velocity 
distribution with a Maxwellian distribution resulted in a 
mean thermal velocity Uth,o = \J'2kT jm = 0.66 mm/s, 
for the unperturbed system. 
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FIG. 2. (color online) Longitudinal (a) and transverse (b) 
current fluctuation spectra from the experiment (color / 
grayscale map) with overlaid Yukawa (k = 1.2) lattice dis- 
persion curves for the two principal lattice directions (lines). 

The determination of the parameters characteriz- 
ing the interparticle interaction (particle charge, Debye 
screening length, plasma frequency) is based on the de- 
pendence of the longitudinal and transverse wave dis- 
persion relations on these quantities, as introduced in 
|24| . Via tracing the trajectories of the particles we cal- 
culate the longitudinal and transverse current fluctua- 
tion spectra from the ri(i) position and Vi(t) velocity 
data as described in [2^. For our system we adopt the 
Debye-Hiickel (Yukawa) type interaction potential, giv- 
ing a pair interaction energy: \J(r) = Q^{r). Match- 
ing the dispersion relations obtained from these spec- 
tra with theoretical dispersion curves of Yukawa lattices 
[26j (see Fig. [5]) we have obtained the following sys- 
tem parameters (within an uncertainty of ±8 %): Debye 
screening length Xd = 0.263 mm, screening parameter 
K = a/Xo = 1.2, particle charge Q — 484Ge (where e 
is the electron charge), and nominal plasma frequency 
loq = \/ nQ"^ jleQTna = 72.1 rad/s. 



B. Experiment with static shear 

As explained in Section II, in this experiment the po- 
larizer P2 (see Fig. [1]) is not rotated, it is set in a way 
that one of the laser beams carries the full power. Two 
methods are available to obtain viscosity data from this 
type of experiments. 

Method 1 is based on the solution of the Navier-Stokes 
(NS) equation. In Q this method was applied for a 
system sheared by two (horizontally shifted), parallel, 
counter-propagating laser beams. We apply the same 
approach for our system sheared by a single beam. In 
this case the NS equation results in a velocity profile 
Vx{y) — Wo exp(— |y|), where vq is the stationary 
flow velocity in the beam center, v is the effective dust- 
background collision frequency (frictional drag), and p 
is the 2D particle mass density. The advantage of this 
method is, that it is relatively noise insensitive and does 
not assume any particular form for the interaction pair 
potential. The drawbacks are that (i) one obtains spa- 
tially averaged data, and (ii) only for the ratio of viscosity 
and collision frequency, thus an estimation or indepen- 
dent measurement of v is needed, which can introduce 
additional uncertainty. 
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FIG. 3. velocity profiles for different shearing laser 

power densities in the dust particle layer in units of the aver- 
age thermal velocity Wth = 1-8 mm/s. Different lines belong 
to the two sides (y > and y < 0) of the sheared region. 

Figure [3] shows time averaged experimental (y) ve- 
locity profiles for different shearing laser powers. K vjr] 
value of 1.64 x 10^'^ kg~^ is derived from these data. The 
average thermal velocity, estimated form the average pe- 
culiar velocities of particles (defined as the velocities rel- 
ative to the flow) within the sheared region, was found 
to lie between vt\i = 1.6 mm/s and 2.0 mm/s for all 
investigated cases using different shearing laser powers, 
about three times higher compared to the unperturbed 
case. This increase is due to the energy absorbed from 
the shearing laser beam. This effect is called shear in- 
duced melting and was studied for dusty plasma crystals 
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in detail in [27|, |28| . As a result of this heating process 
our system is in the strongly coupled liquid phase within 
the sheared region. In the following we use the average 
value vth = 1-8 mm/s for the presentation of our results. 

Method 2 of data acquisition was used so far in molec- 
ular dynamics simulations only [3, [lO, [HI and is based 
on the measurement of the applied shear stress, or, equiv- 
alently the off-diagonal element of the pressure tensor: 



P L 



9 E/ E/ ''"i^y-^ih^ 



(1) 



where A is the area of the region of interest, Vi^x is the 
X component of the peculiar velocity of particle i, rij^y 
is the y component of the distance vector between parti- 
cle i and j, Fij^^ is the x component of the force acting 
on particle i due to the pair interaction with particle j. 
Summation of i is for particles within while summa- 
tion of i runs over all particles interacting with particle 
i. The viscosity is obtained by calculating 



xy 



7 



(2) 



In our experiment the shear rate strongly depends on 
the y-coordinate, due to the exponential velocity pro- 
file shown in Fig. [3] For further evaluation we set up 
a computational grid along the y-axis with a resolution 
of Ay = 0.16 mm. In this way Pxy{y) and 7(2/) can be 
measured by reducing the evaluation volume [A and the 
summation of i in ([T])] to one grid interval at a time. 
Intermediate results obtained with 32 mW/mm^ laser 
power density shear are shown in Fig.Ufa-d), illustrating 
the evaluation process: the observed and time-averaged 
velocity profile (a) is used to compute its derivative, the 
shear rate, which is given in (b) in normalized units, 
7 = {dv X / dy){a/ vt\i) (with uth = 1-8 mm/s). The spatial 
distribution of the off-diagonal element of the pressure 
tensor follows the trend of the shear rate (c) . The viscos- 
ity is the ratio of the data in (c) and (b), and is shown 
in (d) . Substituting the y-coordinatc with the shear rate 
taken from (b) results in the shear rate dependence of 
the viscosity, shown in (e) in normalized units, fj — ri/rjQ, 
where 770 = mnwoa^ — 1.52 x 10^^^ kg/s. 

We have adapted our molecular dynamics (MD) simu- 
lation code [l^, [2l| to the experimental conditions. This 
non-equilibrium MD method is based on the homoge- 
neous shear algorithm [2^, which has frequently been 
used in studies of the viscosity of liquids. The approach 
is based on the Gaussian thermostated SLLOD equations 
of motion and applies the Lees-Edwards sliding periodic 
boundary conditions (for details, see [1^). The shear 
rate, which can be stationary or oscillatory, is an input 
parameter, the off-diagonal element of the pressure ten- 
sor is computed based on ([1]), where the summation vol- 
ume A is the whole simulation box. Viscosity is than 
calculated using ([2]), as it is done for the experiment. 

The simulation results are shown in Fig. Sfe). An im- 
portant difference between the experiment and the sim- 
ulation is in the working principle of the thermostat. In 
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FIG. 4. (color online) Spatial dependence of the (a) velocity 
profile Vxiy); (b) normalized shear rate 7, ; (c) off-diagonal 
element of the pressure tensor —Pxy\ (d) normalized viscosity, 
(e): shear rate dependence of the viscosity, fj{')^^'^), for differ- 
ent laser power densities. The solid lines show MD simulation 
results for F values indicated. 



the simulation a linear velocity profile is achieved to- 
gether with a uniform temperature distribution in the 
whole simulation cell [2^ . In this case the system can be 
well characterized with the Coulomb coupling parameter 
r. In the experiment the temperature is not uniform, 
the system is in crystalline state outside the sheared re- 
gion and is in liquid state inside. An estimation of an 
average (F) w 130 can be given based on the Wth data 
averaged over the investigated region in space. This re- 
sult confirms non- Newtonian behavior [23|, the so called 
shear-thinning effect, where the viscosity decreases with 
increasing shear rate. 

In earlier experiments, like [H-Q, a variation of the 
shear viscosity with changing laser intensity was already 
observed, however, it was attributed to the temperature 
variation near the shearing laser beam(s). The agreement 
between our spatially resolved experimental results and 
our fully thermostated simulations points out, that the 
effect of temperature in this range is marginal, it is the 
high shear rate that causes most of the variations of the 
shear viscosity. 

Taking the average viscosity value of « 0.15 (77 ~ 
2.2 X 10~^^ kg/s) for intermediate shear rates and com- 



5 




■J- 7l/2 71 331/2 271 



° (t) [rad] 

X 



" 0.06 


1 


1 


cogh = 34 rad/s - 


Q.^ 0.04 
0.02 




— 


Y • - 

Pxy ■ 











-0.02 , 








-0.04 








-0.06 


. (b) 

1 


1 


1 



7l/2 71 371/2 271 

(]) [rad] 



FIG. 5. (color online) Measured phase-resolved shear rate (7) 
and off-diagonal element of the pressure tensor {—Pxy) for a 
low (a) and a high (b) frequency periodic shear (symbols) and 
least-square fitted sine functions (lines). 

paring it to the ratio h'/rj = 1.64 x 10^'^ kg~^, as ob- 
tained when applying Method 1, we estimate the ef- 
fective dust-background collision frequency to be z/ « 
3.7 s"^. This value is slightly higher than the ap- 
proximated dust-neutral collision frequency defined as 
i^dn — 3{4:TT/3)r'^nnVnmn/md ~ 3.1 s~^, where and 
nid are the dust radius and mass; n„, Vn and m„ are the 
neutral number density, average speed and mass, respec- 
tively, (5 = 1.26 [13. 

C. Experiment with harmonic shear 

Turning on the rotation of the polar filter P2 (see 
Fig. [1]) results in a sinusoidally modulated periodic shear 
in the x direction. The amplitude of the modulated laser 
power density is set to 32 mW/mm^. The frequency is 
varied between Wsh = 3.7 and 44 rad/s in 25 steps. 

To optimize the signal to noise ratio phase resolved 
averaging is performed: individual oscillation cycles are 
identified based on the reference signal (as explained in 
section II). The period time is divided into ten time slots, 
and every snapshot is assigned to one of them according 
to its actual phase angle cj) relative to the reference sig- 
nal. Despite the low number of particles in narrow slices 
in space and time, this type of averaging ^ 1800 snap- 
shots at a given frequency makes the evaluation Method 
2 (based on the pressure tensor and shear rate) possible. 
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FIG. 6. (color online) Frequency dependence of the shear 
viscosity: (a) magnitude |??(i^)|, (b) complex argument ip{(jj), 
and (c) the real and imaginary parts r;'(a;) and 77" (oj). The 
lines show MD simulation results for F — 200. 

Further, 7(0, y) and Pxy{4>,y) are averaged in space in 
the sheared region 0.5 mm < \y\ < 1.5 mm (see Fig. U]). 
Examples of the measured and Pxy{4') for selected 
frequencies are presented in Fig. [5] together with least 
squares fits in the form /(</)) — ^sin(0 + 0o). 

Performing the sine function least-square fitting pro- 
cedure we obtain the amplitude and phase for both the 
pressure and shear rate for each excitation frequency: 
^■^(wsh), C''(wsh), 0^(wsh), and cj)l{ujsh), respectively. 

The magnitude of the frequency dependent shear vis- 
cosity can be calculated from the amplitudes, while the 
complex argument is given by the difference of the phases: 



as shown in fig. El^a) and (b). 

The real and imaginary (dissipative and elastic) parts 
of the complex viscosity 77(wsh) = '7'('^sh) — i?7"(wsh) are 
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computed using the magnitude and the complex argu- 
ment as 

Vi^sh) = \rii^sh) \ cos [(/^(wsh)] , (4) 
ri" i^sh) = |77(wsh)| sin [(/^(wsh)] • 

With increasing we observe a decrease of rj' and 
an increase of rj" , as shown in Fig. ^c). The crossover 
of the real and imaginary parts can be observed at a 
frequency Wcross/i^o = 0-3 ± 0.05. Our MD simulation 
(using input parameters obtained from the experiment) 
show remarkably good agreement with the experimental 
data. 



IV. SUMMARY 

We have studied the static and dynamic shear viscosity 
of a complex plasma layer via the combination of differ- 
ent experimental and simulation techniques. We have 
developed an optical setup by means of which we applied 
a static, as well as a harmonically modulated shear on 
the particle suspension. 

In the experiments using static shear the viscosity was 
measured via the calculation of the pressure tensor ele- 
ments form the recorded particle trajectories. The mea- 
surements have quantified the shear thinning effect, the 
decrease of the shear viscosity with increasing shear rate. 
These results are consistent with an earlier experiment 
0. Molecular dynamics simulations, modeling a homo- 
geneously sheared liquid with system parameters taken 
from the experiment support our findings and provide a 
possibility to rule out the dominance of thermal effects 
on the observed variations of the shear viscosity. 

In the experiments with sinusoidally modulated shear 
the pressure tensor elements have been measured in a 



time-resolved manner. This way we were able to deter- 
mine the viscoelastic response of the dust layer. The 
data have indicated a phase delay between the pertur- 
bation and the response, that increases with increasing 
frequency: we observed a gradual decrease of the viscous 
term (real part of the viscosity) and a gradual increase of 
the elastic term (imaginary part of the viscosity). These 
experimental results are very similar to those found for 
3D Yukawa liquids in molecular simulations [2l| . 

The appearance of such a rich spectrum of non- 
Newtonian behavior observed in complex plasma experi- 
ments can perhaps be explained by referring to the many- 
fold of interactions and effects between the dust and the 
background plasma. Similarly complex plastic response 
is observed in molecular fluids, like polymer or organic 
solutions, paints, etc. The fact that our molecular dy- 
namics simulations, considering an idealized Yukawa sys- 
tem (i.e. neglecting friction, ion wakes, etc.), well re- 
produce the viscoelastic properties of the experimental 
system, however raises new questions. First, what mi- 
croscopic feature is minimally necessary to result in a 
non-Newtonian macroscopic response? Second, at which 
point does the real complexity of dusty plasmas introduce 
qualitatively new macroscopic features, not reproducible 
by the simplified numerical models based on isotropic 
Yukawa interaction? We can point out only two fea- 
tures, namely, the long range microscopic interaction and 
the strongly coupled (correlated) macroscopic state of the 
model system that may in some way be responsible for 
the computationally observed non-linearities. 
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